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We report the results of a search for the decay B° — > K* + K*~ with a sample of 454 ± 5 million 
BB pairs collected with the BABAR detector at the PEP-II asymmetric-energy e + e~ collider at the 
Stanford Linear Accelerator Center. We obtain an upper limit at the 90% confidence level on the 
branching fraction for B (B° -> K* + K*~) < 2.0 x 10~ 6 , assuming the decay is fully longitudinally 
polarized. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



The study of the branching fractions and angular dis- tributions of B meson decays to hadronic final states 
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without a charm quark probes the dynamics of both 
weak and strong interactions, and plays an important 
role in understanding CP violation. Improved experi- 
mental measurements of these charmless decays, com- 
bined with theoretical developments, can provide sig- 
nificant constraints on the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix parameters [lj and uncover evidence for 
physics beyond the standard model 0, 01 . 

QCD factorization models predict the angular distri- 
bution of the decay of the B meson to two vector par- 
ticles (VV), as measured by the longitudinal polariza- 
tion fraction to be ~ 0.9 for both tree- and penguin- 
dominated decays [J]. Two measurements of the pure 
penguin VV decay B -> cj)K* give f L = 0.52 ±0.08 ±0.03 
and f L = 0.49 ± 0.05 ± 0.03 [5], while f L = 0.81+g;^ ± 
0.06 has recently been measured for the decay B° — > 
K*°K*° [||. Several attempts to understand the values 
of Jl within or beyond the standard model have been 
made Q. Further information about decays related by 
SU(3) symmetry may provide insights into this polariza- 
tion puzzle and test factorization models. 

The decay B° — > K* + K*~ is expected to occur 
through a b — ► u quark transition via VF-exchange, as 
shown in Figure [T] or from final-state interactions. Its 
branching fraction is expected to be small, with Beneke, 
Rohrer and Yang 0] predicting (0.09to:SI£o.io) x 10 ~ 6 ' 
while Cheng and Yang quote (0.1 ± 0.0 ± 0.1) x 10~ 6 , 
both based on QCD factorization. The current exper- 
imental upper limit on the branching fraction at the 
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assum- 



90% confidence level (C.L.) is 141(89) x 10 
ing a fully longitudinally (transversely) polarized system. 
Searches for the related decay B° — » K + K~ have pro- 
duced upper limits on the branching fraction at the 90% 
C.L. in the range (0.4 - 0.8) x 10~ 6 0. 
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FIG. 1: The b -> u IF-exchange diagram for B° -> K* + K*~ 



We report on a search for the decay mode B° — > 
K*+K*-, where K*± refers to the K* ± (892) reso- 
nance, without explicit consideration of interference from 
higher mass K* states, and place an upper limit on the 
branching fraction. Charge-conjugate modes are implied 
throughout and we assume equal production rates of 
B+B- and B°B°. 

This analysis is based on a data sample of 454 ± 5 mil- 
lion BB pairs, corresponding to an integrated luminos- 
ity of 413fb _1 , collected with the BABAR detector at the 



PEP-II asymmetric-energy e + e~ collider operated at the 
Stanford Linear Accelerator Center. The e + e~ center- 
of-mass (cm.) energy is y/s — 10.58 GeV, corresponding 
to the T(4S) resonance mass (on- resonance data). In ad- 
dition, 41.2 fb _1 of data collected at 40 MeV below the 
T(4S) resonance (off-resonance data) are used for back- 
ground studies. 

The BaBAR detector is described in detail in Ref. [13] • 
Charged particles are reconstructed as tracks with a 5- 
layer silicon vertex detector (SVT) and a 40-layer drift 
chamber inside a 1.5-T solenoidal magnet. An electro- 
magnetic calorimeter (EMC) comprising 6580 CsI(Tl) 
crystals is used to identify electrons and photons. A ring- 
imaging Cherenkov detector (DIRC) is used to identify 
charged hadrons and to provide additional electron iden- 
tification information. The average K-tt separation in 
the DIRC varies from 12c at a laboratory momentum of 
1.5GeV/c to 2.5a at 4.5 GeV/c. Muons are identified by 
an instrumented magnetic-flux return (IFR). 

The B° — » K* + K*~ candidates are reconstructed 
through the decay of both K*^ to K®n or with one 
K* ± decaying to and the other to K^n". The 

differential decay rate, after integrating over the angle 
between the decay planes of the vector mesons, for which 
the acceptance is uniform, is 



d 2 T 



1 



r dcos#i<icos( 



sin 2 6>i sin 2 6 2 + f L cos 2 9 1 cos 2 9 2 , ( 1 ) 



where d\ and 6*2 are the helicity angles of the K* + and 
K*~, defined as the angle between the daughter kaon 
(K^ or K^) momentum and the direction opposite to 
the B meson in the K*^ rest frame 



11] 



The charged particles from the K*^ decays are re- 
quired to have at least 12 hits in the drift chamber and a 
transverse momentum greater than 0.1 GeV/c. The par- 
ticles are identified as either charged pions or kaons by 
measurement of the energy loss in the tracking devices, 
the number of photons recorded by the DIRC and the 
corresponding Cherenkov angle. These measurements 
are combined with additional information from the EMC 
and IFR detectors, where appropriate, to reject electrons, 
muons, and protons. 

The K® is reconstructed through its decay to ir + tt~ . 
The candidates are required to have a reconstructed 
mass within 0.01 GeV/c 2 of the nominal K® mass 12], a 
decay vertex separated from the B meson decay vertex by 
at least twenty times the uncertainty in the measurement 
of the vertex position, a flight distance in the transverse 
direction of at least 0.3 cm, and the cosine of the angle 
between the line joining the B and K® decay vertices and 
the K^ momentum greater than 0.999. 

We reconstruct the it through the decay n° — * 77. In 
the laboratory frame, the energy of each photon from 
the 7T° candidate must be greater than 0.04 GeV, the 
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energy of the 7r° must be greater than 0.25 GeV, and 
the reconstructed ir° invariant mass is required to be 
0.12 < m 77 < 0.15GeV/c 2 . 

We require the invariant mass of the K* + candidates to 
be 0.792 < m K7r < 0.992 GeV/c 2 . A B meson candidate 
is formed from the K*^ candidates, with the condition 
that the K*^ candidates originate from the interaction 
region. 

B meson candidates are characterized kinemat- 
ically by the energy difference AE = Eg — y/s/2 
and the beam energy-substituted mass tties = 
[(s/2 + p, -p B ) 2 /Ef -p 2 B ] 1/2 , where {E u p t ) and 
(Eb,Pb) ar e the four-momenta of the T(4S) and B 
meson candidate, respectively, and the asterisk de- 
notes the T(45 l ) rest frame. For a final state with a 
7T°, the total event sample is taken from the region 
-0.1 < AE < 0.2 GeV and 5.25 < m E s < 5.29GeV/c 2 ; 
with no 7T°, the signal AE has a smaller width and the 
region —0.08 < AE < 0.15 GeV is used. The asymmetric 
AE criteria are applied to remove backgrounds from 
charm decays which occur in the negative AE region. In 
both cases, events outside the region \ AE\ < 0.07 GeV 
and 5.27 < toes < 5.29 GeV/ c 2 are used to characterize 
the background. 

We suppress the background from decays to charmed 
states by forming the invariant mass, raj), from com- 
binations of three out of the four daughter particles' 
four-momenta. The event is rejected if 1.845 < ran < 
1.895 GeV/ c 2 and the charge and particle type of the 
tracks are consistent with a decay from a D meson. We 
reduce backgrounds from B° — > (f>K*° by assigning the 
kaon mass to the pion candidate and rejecting the event 
if the combined invariant mass of the two charged tracks 
is between 1.00 and 1.04 GeV/ c 2 . Finally, to reduce the 
continuum background and avoid the region where the re- 
construction efficiency falls off rapidly for low momentum 
tracks, we require the cosine of the helicity angle of the 
K*^ candidates to be in the range — 1.0 < cos(#) < 0.9 
for states without a ir° and —0.9 < cos(0) < 0.9 for de- 
cays with a 7T°. 

To reject the dominant background consisting of light- 
quark qq (g = u,d,s,c) continuum events, we require 
| cos 9t I < 0.8, where 9t is the angle, in the cm. frame, 
between the thrust axis 13( of the B meson and that 
formed from the other tracks and neutral clusters in the 
event. Signal events have a flat distribution in |cos#t|, 
while continuum events peak at 1. 

We use Monte Carlo (MC) simulations of the signal 
decay to estimate the number of signal candidates per 
event. After the application of the selection criteria, the 
average number of signal candidates per event is pre- 
dicted to be 1.08 (1.02) for fully longitudinally (trans- 
versely) polarized decays with no n° in the final state 
and 1.18 (1.10) for decays with one ir° in the final state. 
A single candidate per event is chosen as the one whose 



fitted decay vertex has the smallest \ 2 - MC simulations 
also show that up to 7% (2.4%) of longitudinally (trans- 
versely) polarized signal events with no 7r° are misrecon- 
structed, with one or more tracks originating from the 
other B meson in the event. In the case of signal events 
with one 7r°, the number of misreconstructed candidates 
is 11% (4.3%) for longitudinally (transversely) polarized 
signal events. 

We create a Fisher discriminant T to be used in the 
maximum-likelihood (ML) fit, constructed from a linear 
combination of five variables: the polar angles of the B 
meson momentum vector and the B meson thrust axis 
with respect to the beam axis, the ratio of the second- 
and zeroth-order momentum-weighted Legendre polyno- 
mial moments of the energy flow around the B meson 
thrust axis in the cm. frame [l4j]> the flavor of the other 
B meson as reported by a multivariate tagging algo- 
rithm [15| . and the boost-corrected proper-time differ- 
ence between the decays of the two B mesons divided by 
its variance. The second B meson is formed by creating 
a vertex from the remaining tracks that are consistent 
with originating from the interaction region. The Fisher 
discriminant is trained using MC for signal and qq contin- 
uum MC, off-resonance data and data outside the signal 
region for the background. 

We use an extended unbinned ML fit to extract the sig- 
nal yield and polarization simultaneously for each mode. 
The extended likelihood function is 



£= ^ exp (^" 3 )n 



3 



n j'^'j ( X i'i a j) 



(2) 



We define the likelihood L L for each event candidate i as 
the sum of n,jVj(xi]dj) over three hypotheses j (signal, 
qq background and BB backgrounds as discussed below), 
where Vj (sf< ; dj ) is the product of the probability density 
functions (PDFs) for hypothesis j evaluated for the i-th 
event's measured variables the yield for hypoth- 

esis j, and N is the total number of events in the sample. 
The quantities dj represent parameters in the expected 
distributions of the measured variables for each hypoth- 
esis j. Each discriminating variable a?j in the likelihood 
function is modeled with a PDF, where the parameters 
dj are extracted from MC simulation, off-resonance data, 
or (mESj AE) sideband data. 

The seven variables X4 used in the fit are mESj AE, 
J 7 , and the invariant masses and cosines of the helicity 
angle of the two K candidates. Since the correlations 
among the fitted input variables are found to be on aver- 
age ~ 1%, with a maximum of 5%, we take each Vj to be 
the product of the PDFs for the separate variables. The 
effect of neglecting correlations is evaluated by fitting 
ensembles of simulated experiments in which we embed 
signal and background events randomly extracted from 
fully-simulated MC samples. Any observed fit bias is 
then subtracted from the fitted yield. 
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For the final state with no n°, the two invariant mass 
and helicity angle distributions for each K*^ meson are 
indistinguishable and so we use the same PDF parame- 
ters for both K*^ candidates; for the final state with a 
7T°, we use separate PDFs for K* ± — ► K^ir^ and K* T — > 
K t tt°. For the signal, we use a relativistic Breit-Wigner 
for the K*^- invariant mass and a sum of two Gaussians 
for tties and Ai?. The longitudinal (transverse) helicity 
angle distributions are described with a cos 2 9 (sin 2 0) 
function corrected for changes in efficiency as a function 
of helicity angle. The correction also accounts for the re- 
duction in efficiency at a helicity of ~ 0.78 introduced in- 
directly by the criteria used to veto D mesons. The BB 
backgrounds use an empirical non-parametric function 
for AE, the masses and helicity angles. The continuum 
and the BB background ttjes shapes are described by the 
function xyl — x 2 exp[— f (1 — a; 2 )] (with x — ities/E^ 
and £ a free parameter) [16| and a first- or third-order 
polynomial is used for AE and the helicity angles, re- 
spectively. The continuum invariant mass distributions 
contain real K candidates; we model the peaking mass 
component using the parameters extracted from the fit 
to the signal invariant mass distributions together with 
a second-order polynomial to represent the non-peaking 
component. The Fisher distributions are modeled using 
an asymmetric Gaussian for all hypotheses. 

BB backgrounds that remain after the event selection 
criteria have been applied are identified and modeled us- 
ing MC simulation based on the full physics and detector 
models [TtJ • There are no significant charmless BB back- 
grounds. The charm BB backgrounds are effectively sup- 
pressed by applying the veto on D meson mass described 
above. The remaining charm BB background events are 
mostly single candidates formed from the decay products 
of a D, D* or D^, together with another track from 
the event. Given the uncertainty in the polarization and 
branching fractions of these backgrounds, we allow the 
BB background yield to float in the fit. 

The continuum background PDF parameters that are 
allowed to vary are the T peak position, £ for jtieSj the 
slope of Ai?, and the polynomial coefficients and normal- 
izations describing the mass and helicity angle distribu- 
tions. We fit for the branching fraction B and fi directly 
and exploit the fact that B is less correlated with fjj than 
is either the yield or efficiency taken separately. We val- 
idate the fitting procedure and extract fitting biases by 
applying the fit to ensembles of simulated experiments 
using the extracted fitted yields from data. The qq com- 
ponent is drawn from the PDF, and the signal and BB 
background events are randomly sampled from the fully 
simulated MC samples. 

The total event sample consists of 602 and 1923 events 
for B° — > K* + K*~ with zero or one ir° in the final 
state, respectively. The corresponding signal event yield 
is 1.8l" 2 '7 an d ^.lig'g ano - ^ ne longitudinal polarization 
jh is 0.0 ±0.6 and 1.0 ±1.0, respectively. Given the 



large errors on fi,, we repeat the analysis with set 
to 1.0; this gives the most conservative 90% confidence 
level upper limit on the branching fractions. The results 
of the ML fits with Jl = 1.0 are summarized in Table Q] 
The BB background yield agrees with the MC predic- 
tion within the large statistical errors. We compute the 
branching fractions B by subtracting the ML fit bias from 
the fitted yield and dividing the result by the number of 
BB pairs and by the reconstruction efficiency, e, times 
B(K° -> K° -> 7T+7T-) = 0.5 x (69.20 ± 0.05)% and 
B(tt° -> 77) = (98.80 ± 0.03)%. The significance S of 
the signal is defined as S = 2Aln£, where Aln£ is the 
change in likelihood from the maximum value when the 
number of signal events is set to zero, corrected for the 
systematic errors defined below. 

The significance of the B° — ► K* + K*~ branching frac- 
tion is 0.87cr, including statistical and systematic uncer- 
tainties. The 90% C.L. branching fraction upper limit 
(2?ul) is determined by combining the likelihoods from 
the two fits and integrating the total likelihood distri- 
bution (taking into account correlated and uncorrelated 
systematic uncertainties) as a function of the branching 
fraction from to B UL , so that f,f UL CdB = 0.9 / °° CdB. 

Figures [2] and [3] show the projections of the two fits 
onto mES, AE, K*^ mass and cosine of the K*^ helicity 
angle for the final state with zero and one ir° , respectively. 
The candidates in the figures are signal-enhanced with a 
requirement on the probability ratio V a i g /(V s i s + ^bkg), 
optimized to enhance the visibility of potential signal, 
where "P S i g and Vbkg are the signal and the total back- 
ground probabilities, respectively (computed without us- 
ing the variable plotted). The dip in helicity at ~ 0.78 is 
created by the criteria used to veto charm background. 

TABLE I: Summary of results with /l=1.0 for the fitted 
yields, fit biases, reconstruction efficiencies e, sub-branching 
fractions 11^1 branching fraction B (B° — > K* + K*~), sig- 
nificance S, and 90% C.L. upper limit Bul- The first error is 
statistical and the second, if given, is systematic. 



Final State 


K° n+ K° s TT" 


K° TT* TT 


Yields (events): 






Total 


602 


1923 


Signal 






BB bkg. 


20 ±20 


84 ±51 


qq bkg. 


580 ± 25 


1835 ± 70 


ML Fit Biases 


-0.170 


1.70 


Efficiencies and B: 






e(%) 


8.89 ±0.08 


4.83 ±0.04 




5.32 


15.19 


B (xl0~ 6 ) 


0.38+J J ±0.05 


0.76+^ ±0.16 


Significance S (a) 


0.50 


0.74 


Combined Results: 
B (xbT 6 ) 


fi r 9 +0.S3+0.08 
u - dz -0.58-0.06 



\ / — U.vJO" 

Significance S (a) 0.87 



g UL (xl0~ 6 ) 




0.9 ~ 0M 
KMassfGeV/c 2 ) 

FIG. 2: Projections of the multidimensional fit onto (a) to.es; 
(b) AE; (c) K mass; and (d) cosine of K helicity an- 
gle for B° -> K m± (-> Kl-K^K*^ '(-* K%ir*) events selected 
with a requirement on the signal-to-total likelihood probabil- 
ity ratio, optimized for each variable, with the plotted variable 
excluded. The points with error bars show the data; the solid 
line shows signal-plus-background; the dashed line is the con- 
tinuum background; the hatched region is the signal; and the 
shaded region is the BB background. 




0.9 0.95 

K ^K'V MassfGeV/c 2 ) K"->IC'V helicity 

FIG. 3: Projections of the multidimensional fit onto (a) toes; 
(b) AE; (c) K mass; and (d) cosine of K* ± helicity angle 



for B 



i) 



K 



<±t 



K*^(~^ K^ty ). The same projec- 



tion criteria and legend are used as in Fig. [2] 



The systematic uncertainties are summarized in Ta- 
ble lU The errors on the branching fractions arise from 
the PDFs, fit biases and efficiencies. The PDF uncer- 
tainties are calculated by varying the PDF parameters 
that are held fixed in the original fit by their errors. The 
uncertainty from the fit bias includes its statistical un- 
certainty from the simulated experiments and half of the 
correction itself, added in quadrature. The uncertain- 
ties in PDF modeling and fit bias are additive in nature 
and affect the significance of the branching fraction re- 
sults. Multiplicative uncertainties include reconstruction 
efficiency uncertainties from tracking and particle iden- 
tification (PID), track multiplicity, MC signal efficiency 



statistics, and the number of BB pairs. 



TABLE II: Estimated systematic errors in the final fit. Error 
sources which are correlated and uncorrelated when combined 
from the two decays are denoted by C and U, respectively. 



Final State 


K° s n+ K° s TT" 


1^0 TS± „° 
A s IT A 7T 


Additive errors (events) 






Fit Bias [U] 


0.09 


0.85 


lit Parameters [U] 


0.06 


0.25 


Total Additive (events) 


0.10 


0.88 


Multiplicative errors (%) 






Track Multiplicity [C] 


1.0 


1.0 


MC Statistics [U] 


0.5 


0.6 


Number of BB pairs [C] 


1.1 


1.1 


PID [C] 


2.2 


1.1 


Neutrals Corrections [C] 




3.0 


K° s Corrections [C] 


1.8 


1.4 


Tracking Corrections [C] 


1.6 


0.8 


Total Multiplicative (%) 


3.6 


3.9 


Total B error (xl(T b ) 


0.05 


0.16 



In summary, we have measured the branching frac- 



tion B{B° -> K*+K 



[0.52 



+0.83+0.08] 
-0.58-0.06J 
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as- 

suming the decay is fully longitudinally polarized. The 
90% C.L. upper limit on the branching fraction B(B° — ■> 
K* + K*~) < 2.0 x 10~ 6 is nearly two orders of magnitude 
more stringent than previous searches. 
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